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Measurement-based quantum computation has revolutionized quantum information processing, 
and the physical systems with which it can be implemented. One simply needs the ability to 
prepare a particular state, known as the cluster state, and subsequently to perform single-qubit 
measurements on it. Nevertheless, a scalable implementation is yet to be realized. Here we propose 
a hybrid light-matter system comprised of coupled cavities interacting with two level systems. Util- 
ising the stable, individually addressable, qubits resulting from the localised long-lived atom-photon 
excitations, we demonstrate how to use the natural system dynamics to 'weave' these qubits into 
a cluster state and propose the implementation of quantum algorithms employing just two rows 
of qubits. Finally, we briefly discuss the prospects for experimental implementation using atoms, 
quantum dots or Cooper pair boxes. 



Quantum computation on cluster states [l[ has been 
proposed in a variety of systems, including linear optics, 
quantum dots, neutral atoms in optical lattices, and fly- 
ing atom schemes @, S i, i, @, 0, 1, H To date, 
experiments have been performed using optical lattices 
where the cluster state can be created, but the cur- 
rent lack of individual addressing remains the stumbling 
block and linear optics 13, IH , where scalability remains 
a problem due to the need to generate the initial many- 
photon state from, for example, high orders of the para- 
metric down conversion proecss. On the other hand, 
there have recently been theoretical and experimental 
breakthroughs into the possibility of direct coupling of 
high Q cavities and in achieving strong coupling between 
the cavity mode and an embedded two-level system. A 
variety of technologies have been employed, namely fiber 
coupled micro-toroidal cavities interacting with atoms 
[3, EH , arrays of defects in photonic band gap materi- 
als (PBGs) 3, 17, 18 1 and superconducting qub its cou- 
pled through microwave stripline resonators [l9(. This 
has prompted proposals for the implementation of opti- 
cal quantum computing [2(| , the production of entangled 
photons 2l| and the realization of Mott insulating and 
superfluid phases 2^, HI 24| . Here we propose the use of 
such arrays for the realization of cluster state quantum 
computation. 

System Description: We start be describing the sys- 
tem and showing how to construct qubits from the hy- 
brid light-matter excitations (polaritons). For simplicity, 
we describe the system as a linear chain of N coupled 
cavities doped with two level systems, although this is 
readily adapted to the two-dimensional setting that we 
require. | g) k and | e) k are the atomic ground and excited 
states at site k (we henceforth use the term 'atom' to re- 
fer to any relevant two level system). The Hamiltonian 
describing the system is the sum of three terms; H* ree 
is the Hamiltonian for the free light and dopant parts, 
H mt the Hamiltonian describing the internal coupling of 
the photon and dopant in a specific cavity and H hop for 



the light hopping between cavities. 
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Ud and A are the photon frequencies and hopping 
rates respectively and g is the light-atom coupling 
strength. The H' ree + H mt component of the Hamil- 
tonian can be diagonalized in a basis of combined pho- 
tonic and atomic excitations, called polaritons (Fig. [T]). 
These polaritons are defined by creation operators 
P k ' n ^ = | n±) k (g, | fc , where the polaritons of the A:th 
atom-cavity system are given by |n±) fe = (\g, n) k ± 
| e, n — l) fc )/y2 with energies E^ = nuidig^/n, and | n) k 
denotes the n-photon Fock state. As has been shown 
elsewhere, a polaritonic Mott phase exists in this system 
where a maximum of one excitation per site is allowed 
[2^ | . This originates from the repulsion due to the photon 
blockade effect 2|| 2||. In this Mott phase, the system's 
Hamiltonian can be written in the interaction picture as 
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fc x k+1 + P k P ] k+1 , where P fe T = x k 
(Fig. [TJ . As double or more occupancy of the sites is pro- 
hibited, one can identify P^ with a k = a k +ia k , where a k 
and a\ are the standard Pauli operators. The system's 
Hamiltonian then becomes the standard XY model of in- 
teracting spin qubits with spin up/down corresponding 
to the presence/absence of a polariton. 
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Some applications of XY spin chains in quantum informa- 
tion processing can thus been implemented in this system 
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FIG. 1: We work with a 2D array of atom-cavity systems. When the atom is on resonance with the cavity, the ground state 
\g, 0) and the first excited state |1— ) of the combined atom-photon (polaritonic) system in each site can be used as qubits. By 
applying Stark shifts with control electrodes or properly tuned laser fields to sets of qubits (the gates shown under the qubits), 
we disable the exchange Hamiltonian of a qubit to all of its neighbours. Applying gates B, C, D or equally switching off gate 
A(part (a)) creates chains of 3 qubits which apply controlled-phases and SWAPs between the qubits at either end of the chain 
(indicated by dashed lines) . Repeating the sequence with the other 3 gates is sufficient to connect the 3-chains and generate a 
cluster state in parallel across the whole device. Single qubit rotations and measurements are made by properly applying local 
external fields, utilizing the fact that the cavities can be well separated. 



Cluster state generation: The typical implementation 
of cluster state quantum computing requires initializing 
all qubits in a 2D lattice in the | +) = (| 0) + | 1))/V2 
state and then performing controlled-phase gates (CP) 
between all nearest-neighbours. In the present system, 
we have no direct two-qubit gate and the available inter- 
action is not of the Ising type, which straightforwardly 
gives controlled-phase gates but an 'always on' global 
Hamiltonian coupling of the XY form. It behooves us to 
consider how we will proceed with the measurement se- 
quence once the cluster state has been generated without 
the system continuing to evolve. (The need to disable the 
evolution is an aspect often neglected when discussing the 
generation of cluster states, either as the ground states 
of 'natural' Hamiltonians [28|, or from evolution of the 
Ising Hamiltonian) . This requirement can be realized by 
combining the system's natural dynamics with a protocol 
where some of the available physical qubits are used as 
gate "mediators" and the rest as the logical qubits. The 
mediator atoms can be Stark shifted on and off resonance 
from their cavities through the application of an external 
field, inhibiting the photon hopping and thereby isolat- 
ing each logical qubit. The same inhibition of couplings 
will be used to generate the cluster state. We note here 
that the error introduced in the step is due to a second- 
order transition between on-resonance qubits (via a dark- 
passage through the central off-resonant qubit), which is 
thus supressed by a factor of order A/ A, where A is the 
detuning of the off-resonant cavity. 

Before describing the 4-step global gate sequence to 
create the cluster state, first observe that for the control 
phase part is enough to localize chains of 3 qubits, let 
them evolve for a time to = n/(2y/2A) and then apply 
a measurement on the middle 'mediator' qubit (in the 



a z basis). Depending on the measurement result, | 0) or 
| 1), a nonlocal gate is generated between the remaining 
two qubits, either SWAP.(ct z <g> a z ).CP or SWAP. CP 
respectively [29l I30I] . In both cases, the gates in addition 
to the CP are Clifford operations which can be recorded 
and taken into account during the measurement-based 
computation. 

Our sequence to generate the cluster state initiates by 
preparing all qubits in the | +) state through the appli- 
cation of global 7r/2 pulse. One quarter of the sites will 
be used as logical qubits and the rest as "mediators" and 
"off' qubits interchangeably. All qubits addressed by the 
gates A to D (Fig.[TJ) are, by default, "off" , thereby isolat- 
ing all the qubits. Switching on any one of the four gates 
thus creates chains of 3 qubits, which we use to enact 
a CP between pairs of qubits (separated by a mediator 
qubit, which was previously off). Consecutive use of each 
of the gates A to D serves to enact a CP gate between 
a particular qubit and all of its nearest-neighbours, and 
this happens in parallel across the whole device. This 
entire sequence is illustrated in Supplementary Video 1. 
The measurement sequence is then applied as requested 
by the cluster state algorithm, utilizing the local acces- 
sibility of the sites (in any implementation, the cavity- 
atom systems are well separated compared to the res- 
olution of the external field used for addressing them) 
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In Fig. [21 we calculate the fidelity of generation of a 
cluster state on a 3x3 array of cavities. More sophisti- 
cated schemes have the potential to further reduce the 
experimental errors. For example, standard Hamiltonian 
simulation techniques allow us to negate the second or- 
der exchange term due to the off-resonance cavities, sim- 
ply by repeatedly applying o z gates to every second on- 




FIG. 2: The fidelity of generation of a cluster state on a 
3x3 grid of cavities, as the detuning A of the mediator off- 
resonance cavities is varied(in units of the hopping A). The 
dashed line includes post-selection on getting | 0} outcomes 
when measuring off-resonance qubits. The gray lines also 
incorporate spontaneous decay and cavity leakage of 0.05A 
(dark) and 0.08,4 (light). 

resonance triplet throughout the evolution. One might 
even hope that we could use this coherent effect to en- 
hance the scheme through the use of, for example, op- 
timal control techniques. Most of the errors considered 
here (cavity leakage, spontaneous emission of the atom, 
and on-off detuning of qubits) are local effects, introduc- 
ing local noise, which can ultimately be addressed by 
fault-tolerant techniques [3lj ]. 

Implementing algorithms: Initial experimental algo- 
rithmic implementations with coupled cavities can be ex- 
pected to utilize the most basic building block of our 
scheme, a 3 x 3 grid of cavities, which allows us to gen- 
erate a four-qubit cluster state. As with the four-photon 
cluster state recently used by Walther et al. |12j |. this 
cluster state would be suitable for demonstrating the 
preparation of an arbitrary one-qubit state, an entan- 
gling gate between two qubits, and even the implementa- 
tion of Grover's search algorithm on two qubits 12j. For 
example, by applying the local gates H (g> H (g) a z % a z , 
where H is the Hadamard rotation, we convert of 'box' 
cluster that the 3x3 grid prepares into the ID cluster 
state of 4 qubits, which is given the interpretation of a 
single qubit, and measurements on the state yield quan- 
tum gates on this single qubit. Moreover, generation of 
this four qubit cluster state is simpler than generation 
of an arbitrarily sized cluster state because we only need 
two control steps instead of four, thereby keeping us even 
further within the decoherence time of the system. 

Perhaps the next important step would then be to 
demonstrate Shor's factoring algorithm, the factoring of 
15 being the standard demonstration. To implement as a 
cluster state computation, the six computational qubits 
|33| translate into the requirement of a cluster state that 
is eleven qubits wide. Hence, we need an array which is 
21 cavities wide. The breadth of the cluster state, which 
corresponds to time in the circuit model, is a quantity 
that we can trade against the time taken for the compu- 
tation. At one extreme, we can create the whole cluster 



FIG. 3: Sequence for minimising the number of qubits re- 
quired for a cluster state computation, (a) After the first 
n — 1 steps of the algorithm, the first column of qubits is ini- 
tialised in the | +) state, and the third column, with qubits 
denoted by *, are in the state of output for the first n— 1 steps 
of the computation, (b) We use control sequences, bringing 
mediator qubits on resonance, to convert the | +) states into 
a cluster state, and to entangle them with the output qubits. 
The SWAP in the entangling operation moves these output 
qubits to the first column, (c) Measure the qubits of the first 
column as corresponds to the n th step of the computation, 
and reinitialise in the | +) state. The rightmost column cor- 
responds to the output. The sequence then repeats. 

state in one go, with the simple set of four steps already 
outlined, and we benefit from the large degree of paral- 
lelism available to us. This requires a 2D grid of cavities 
of size 21 x 311 34]. At the other extreme, a grid of 21 x 3 
cavities suffices. In this case, one starts with the 11 x 2 
cluster state, and performs one time step of measurement 
(i.e. measure the 11 qubits in one column). The result 
remains in the other column. We then repeat the clus- 
ter state generation process, reinitialising the measured 
qubits in the cluster state, and performing the next time 
step (Fig. [3]). This requires 156 consecutive entangling 
steps, but the reinitialising of the cluster state after mea- 
surement eliminates the effect of decoherence over this 
timescale. Any combination between these two extremes 
is also possible, and is a necessary property of any scal- 
able implementation of cluster state computation for the 
sake of preventing decoherence. 

Once initial cluster state experiments have been per- 
formed, it simply becomes a question of how many cav- 
ities one can reasonably couple together. Alternatively, 
since the two-qubit gate that we can generate is entan- 
gling (and hence universal for quantum computation), we 
can also consider using it directly to implement the cir- 
cuit model of computation. This has a much smaller over- 
head of qubits, but instead requires much higher quality 
cavities. For example, to factor 15 we would only need a 
5x3 grid of cavities to give us six computational qubits. 
However, we would need approximately 15 consecutive 
entangling steps (we have attempted to minimise this 
number by allowing as many of the gates to be applied 
in parallel as possible, and by optimising the initial la- 
belling of each qubit), hence requiring a time of order 
15t:/(V2A). Hence, to reduce the effect of dissipativc 
decay, we require an order of magnitude improvement in 
the decoherence properties of the qubits to compensate 



for the increased running time. 

Experimental implementations: As previously men- 
tioned, there are three primary candidate technologies; 
fiber coupled micro-toroidal cavities 14, [ll|, arrays of 
defects in PBGs [H, [l7l HH and superconductin g qu bits 
coupled through microwave stripline resonators [19(. In 
order to achieve the required limit of no more than one 
excitation per site 22], the ratio between the internal 
atom-photon coupling and the hopping of photons down 
the chain should be of the order of g /A ~ 10 2 — 10 1 (A can 
be tuned while fabricating the array by adjusting the dis- 
tance between the cavities and g depends on the type of 
the dopant). In addition, the cavity/atomic frequencies 
to internal coupling ratio should be cud, ^0 ~ 10 4 g, 10 5 g 
and the losses should also be small, gj maxfre, 7) ~ 10 3 , 
where k and 7 are cavity and atom/other qubit decay 
rates. The polaritonic states under consideration are es- 
sentially unaffected by decay for a time 1Q/A (10ns for 
the toroidal case and 100ns for microwave stripline res- 
onators). While the decay time of 10/vl may seem un- 
comfortably close to the preparation time for a cluster 
state, y/2n/A, the previously described technique (Fig. [3]) 
of continuously reforming the cluster state and connect- 
ing it to the output of the previous stage allows a con- 
tinuous computation that exceeds the decay time for an 
individual cavity. The required parameter values are cur- 
rently on the verge of being realised in both toroidal mi- 
crocavity systems with atoms and stripline microwave 
resonators coupled to superconducting qubits, but fur- 
ther progress is needed. Arrays of defects in PBGs re- 
main one or two orders of magnitude away, but recent 
developments, and the integrability of these devices with 
optoelectronics, make this technology very promising as 
well. In all implementations the cavity systems are well 
separated by many times the corresponding wavelength 
of any local field that needs to be applied in the system 
for the measurement process. 

Conclusions: In this paper, we have shown how 
universal quantum computation could be realized in a 
coupled array of individually addressable atom-cavity 
systems, where the qubits are given by mixed light- 
matter excitations in each cavity site. While single- 
qubit operations can be locally achieved, the only avail- 
able interaction between qubits is due to the natural 
system Hamiltonian. We show how to manipulate this 
to give a controllcd-phase gate between pairs of qubits. 
This allows computation either using the circuit model, 
or a measurement-based computation, the latter being 
most suited to reducing experimental errors. We have 
discussed possible architectures for implementing these 
ideas using photonic crystals, toroidal microcavities and 
superconducting qubits and point out their feasibility 
and scalability with current or near-future technology. 
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a nearest-neighbor, 2-qubit gate algorithm. Hence, the number of qubits remains, 

possibility for some small degree of optimisation in the 



